Finite-time effects and ultraweak ergodicity breaking in superdiffusive dynamics 
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We study the ergodic properties of superdiffusive, spatiotemporally coupled Levy walk processes. 
For trajectories of finite duration, we reveal a distinct scatter of the scaling exponents of the time 
averaged mean squared displacement 8 2 around the ensemble value 3 — a(l<a<2) ranging from 
ballistic motion to subdiffusion, in strong contrast to the behavior of subdiffusive processes. In 
addition we find a significant dependence of the trajectory-to-trajectory average of S 2 as function of 
the finite measurement time. This so-called finite-time amplitude depression and the scatter of the 
scaling exponent is vital in the quantitative evaluation of superdiffusive processes. Comparing the 
long time average of the second moment with the ensemble mean squared displacement, these only 
differ by a constant factor, an ultraweak ergodicity breaking. 

PACS numbers: 87.10.Mn,89.75.Da,87.23.Ge,05.40.-a 



Suppose you are recording the trajectories of individ- 
ual blue sharks in the ocean over time t. Calculating the 
time averaged mean squared displacement (MSD) S 2 for 
each shark you find that some animals move almost bal- 
listically while others appear to move much slower. Does 
this indicate that the animals follow different generic mo- 
tion patterns? As we show here for the celebrated Levy 
walk (LW) model of superdiffusion, the intrinsic non- 
ergodicity in trajectories of finite length indeed gives rise 
to a wide distribution of apparent scaling exponents, even 
to subdiffusive values, although the motion is produced 
from identical distributions. This surprising finding is ac- 
companied by a strong reduction of the amplitude of S 2 
at finite measurement times and strongly contrasts the 
non-ergodicity observed in subdiffusive motion. 

Blue sharks are indeed just one example of marine 
predators followed over large distances that show scaling 
laws in their foraging behavior consistent with LW dy- 
namics [H, similar to findings from other tracking stud- 
ies of individual animals or humans LWs are more 
widely applied, inter alia to describe intermittent chaotic 
systems p-0| , turbulent flow [§] , accelerated diffusion in 
Josephson junctions [j|, negative Hall-resistance in semi- 
conductors [13], diffusion of atoms in optical lattices 11] 
and of light in disordered media [13] , blinking statistics of 
quantum dots 13], movement strategies in mussels 14], 
or even T-cell motility in the brain [l5j. Many of these 
systems may be analyzed on the single trajectory level. 

Despite this ubiquity of LWs their ergodic behavior has 
not been studied in detail. However the question whether 
a system is ergodic becomes relevant when instead of the 
conventional MSD (x 2 (t)) = J x 2 P(x,t)dx defined as en- 
semble average over the probability density P(x, t) we 
use time averages over single trajectories. For time series 
x(t) of duration T the time averaged MSD is defined via 



where r denotes the lag time. The behavior of Sx 2 has 
been studied in detail for the subdiffusion case, (x 2 (t)) ~ 
t 1 with < 7 < 1, revealing distinct discrepancies be- 
tween ensemble and time averaged MSD for scale free 
waiting time processes 0, 17]. This so-called weak er- 
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godicity breaking (WEB) means that (x 2 (t)) ^ Sx 2 even 
for long T [lllH. 

while other subdiffusive processes such 
as fractional Brownian motion are ergodic in the sense 
that (x 2 {t)) = Sx 2 for sufficiently long T [3, [13]. WEB 
has indeed been observed in experiments, for instance, 
for the motion of protein channels in the walls of living 
cells [13] and of lipid granules in yeast cells [2l| . 

To study the ergodic properties of of LWs we re- 
call their definition within continuous time random walk 
(CTRW) theory [13]. A CTRW is based on the joint dis- 
tribution ty(x,t). For each jump we draw from ^(x, t) 
a random waiting time t and jump length x [2314251 ]. 
To describe superdiffusive processes (x 2 (t)) ~ i 7 with 
7 > 1, LWs are endowed with a spatiotemporal cou- 
pling for which we chose the simplest form ^(x, t) = 
^ip(t)S(\x\ — vt) [24j. Confined by an expanding horizon 
at positions ±t>i from the origin, this CTRW performs 
statistically independent free paths with constant veloc- 
ity \v\, whose durations are distributed according to the 
power law ip(t) = J ^(x,t)dx ~ For < a < 1 

the resulting motion is ballistic, 7 = 2, for 1 < a < 2 we 
observe sub-ballistic superdiffusion with 7 = 3 — a, while 
for a > 2 the motion is normal diffusive, 7 = 1 [3, 13] ■ 
The mean sojourn time (t) = J Q tip(t)dt is infinite for 
< a < 1 and finite otherwise. In contrast to Levy 
flights with their diverging variance (26|, LWs are thus 
physical models for particles with a maximum propaga- 
tion speed. Apart from the description in terms of above 
continuous time random walk scheme with fy(x,t), LWs 
can be described as a renewal process [3], in terms of a 
master equation [53] , a fractional transport equation [53] , 
or a Langevin approach based on subordination [29| . 
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Here we focus on the behavior of time averages and 
ergodic properties of LWs in the relevant superdiffusive 
range 1 < a < 2. From analytical results and extensive 
numerical simulations we highlight the particular role of 
the finiteness of trajectories when calculating the time 
averages. Namely, we show that the scaling exponents 
of S 2 apparently become random quantities and that the 
amplitude of the time averages is a function of the mea- 
surement time T. Moreover, we report an ultraweak er- 
godicity breaking of superdiffusive LWs. These effects are 
important to interpret time averages of LW processes. 

A full analytical solution for the the time averaged 
MSD Sx 2 is obtained from the renewal framework (9(. 
The starting point is the velocity autocorrelation func- 
tion C v (t) = limt/^oo \t' - t|— 1 flf^ 1 v(t")v(t" + t)dt", 
where the time average is taken over a trajectory of infi- 
nite length. In the velocity model for LWs employed here 
the velocity fluctuates between +v and — v with equal 
probability, meaning that only single events contribute 
to C v (t), which in turn is the result of an averaging of 
event durations along a trajectory. The problem can be 
rephrased in terms of the probability that a walker is in 
an ongoing event of duration T between and t given 
that we pick an arbitrary origin on the time axis. To 
obtain C v we simply average over all such possible du- 
rations. Once C v is known, 5x 2 is readily obtained from 



the Green-Kubo formula |3f| 5x 2 (t) = 2 /J" '( T -t)C v (t)dt. 
For infinite trajectories, we obtain the result 
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where we have set |u| = 1. As for 1 < a < 2 the 
mean waiting time (t) is finite, individual trajectories at 
sufficiently long (infinite) times become self-averaging, 
such that there will be no difference between Sx 2 (t) ob- 
tained from different trajectories and the trajectory-to- 
trajectory averaged quantity (Sx 2 (t)). In other words, 
the actual series of events is irrelevant in the case of in- 
finite trajectories. In reality one never deals with in- 
finite trajectories, albeit they might become extremely 
long. Once a trajectory is finite, irrespective of its actual 
length, there always exits a non-zero probability that the 
walker will be 'locked' in a single motion event persisting 
along a great fraction or even during the entire trajectory. 
Not surprisingly, the MSDs Sx 2 (t) of individual trajec- 
tories will not coincide but show a scatter of amplitudes, 
as shown below. 

In our simulations we use the concrete form ip(t) = 
ail + for the waiting time distribution. From 

this asymptotic power-law we generate M = 10 4 time se- 
ries of particle coordinates xj(t), where j labels different 
trajectories. We calculate the ensemble averaged MSD 
(Ax 2 (r)} = M _1 ^^ x Xj(r) 2 and the time averaged 
MSD through Eq. (JTj . Fig. [TJ shows typical results for 
Sx 2 for 400 different trajectories of duration T = 10 8 time 




FIG. 1: Time averaged mean squared displacement Sx 2 ob- 
tained from 400 trajectories for a = 1.2 (top) and a — 1.5 
(bottom). Full red lines depict different scaling behaviors as 
indicated, the full black lines are averages over all trajectories. 



steps (A) for a — 1.2 and a = 1.5. Remarkably, while 
Sx 2 for all trajectories coincides and shows superdiffusive 
scaling at shorter lag times, at longer r, Sx 2 displays a 
wide spread of slopes ranging from ballistic motion to 
subdiffusion (7 < 1). At the same time the ensemble- 
averaged MSD predicts a unique long-time scaling of the 
form (x 2 (t)) ~ < 3 ~ Q , confirmed by our simulations (not 
shown here). Thus ergodicity, the equivalence of long 
time and ensemble average is broken. Moreover, self- 
averaging does not take place. In contrast to subdiffu- 
sive CTRW with diverging mean waiting time, where the 
scaling is identical for all trajectories but the generalized 

mini, 



diffusion coefficient becomes a random variable 
here we observe that the scaling exponent of individual 
trajectories appears random. We note that this effect is 
not due to bad statistics at larger t as t C T is fulfilled 
for all t shown in Fig. [TJ Performing an average over 
all trajectories, (Sx 2 ), the full black lines in Fig. [TJ the 
result seems to follow the scaling predicted by Eq. ©. 
However, this agreement is only apparent, see below. 

We first quantify the deviations between different tra- 
jectories in terms of the distribution of S 2 around the 
trajectory-to-trajectory average (5x 2 ), 



5x 2 (t) 

WW) 



(3) 



The results for a — 1.2 and a — 1.5 are shown in Fig. [5] 
(Left) . In the case of an infinite trajectory we would find 
a sharp peak at £ = 1, which is approximately observed 
for the shortest r. In contrast for finite-time trajecto- 
ries P(£|t) apparently relaxes towards a skewed limit- 
ing distribution with a maximum well below the ergodic 
value £ = 1. Therefore the average value appears to be 
dominated by one or few very long waiting time events 
locked onto a given velocity mode, and the self-averaging 
is not fulfilled. In addition we measure the long-time 
scaling in individual trajectories by least-squares fit of 
the last decade of Sx 2 to the power-law t a * vp , obtaining 
the scatter distribution of the apparent scaling exponent 
a app- The resulting distributions -P(a a pp) for a = 1.2 
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FIG. 2: Left: Scatter distribution P(£) of the time averaged 

MSD Sx 2 versus £ = Sx 2 /(8x 2 ) for a = 1.2 (black) and a = FIG. 3: (Sx 2 ) (dashed red line) and (5x 2 ) f (full lines) with 

1.5 (red). Right: Scatter distribution P(a app ) of the apparent t = 10 7 A (black symbols) and T = 10 8 A (blue symbols), 

scaling exponent a app obtained by least-squares fit of Sx 2 ~ The right panel shows the same plots on the logarithmic scale. 
£ aapp to the last decade of the time series in Fig. [1] 



and a = 1.5 are shown in Fig. [2] (Right). We see that the 
maximum of -P(a app ) is well below the infinite-time av- 
erage cv ap p = a, and that a non-negligible fraction of tra- 
jectories in fact exhibits subdiffusion. This demonstrates 
that the time average of a super-diffusive dynamical pro- 
cess can in fact display subdiffusive behavior on the level 
of single trajectories of finite duration. Again, we see 
that finite time averages such as 5 2 are obviously dom- 
inated by either extremely long motion events pushing 
<*app to values closer to a app = 2, while strong oscilla- 
tions between velocity modes induce localization effects 
and values a app < 1. These observations will be crucial 
for the correct interpretation of single trajectory mea- 
surements of superdiffusive processes. 

Having established that there is no unique scaling of 
5x 2 along finite-time single trajectories one might won- 
der whether and how the finiteness of single trajecto- 
ries affects the corresponding average over an ensemble 
of trajectories. This problem can be treated exactly 
with the renewal approach. Once an arbitrary origin 
is specified on the time axis the probability that the 
walker is in a motion event of duration $ at time is 
po(§)dd — $ip('&)d'&/'& ) where i? is the average time span 
of •& along a finite-time trajectory of duration T and 
ensures the correct normalization, $ = J Q T = 
^-(1 - a(l + aT)(l + T)- a ). The probability that the 
event persists until t is p p (£|i?) = (# — such that 

the probability that the walker is in a motion event of 
duration d between and t is po(-d)p p (t\'d)d'd. The ve- 
locity autocorrelation function C[ (t) is then obtained by 
averaging over all possible durations up to T, such that 



for t < T, such that for finite-time trajectories we find 



cUt) 
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For long T the time averaged MSD has the form 



(5) 



Sx 2 (t))+T 



->S— a 



a - 1 



- (6) 

Simulations results for {Sx 2 ) are shown in Fig.[3j demon- 
strating good agreement with the result ((3}. Indeed 
we find that on a logarithmic scale (Right, the conven- 
tional representation of time averaged MSD data) one 
hardly observes deviations from Eq. @, however, on a 
linear scale pronounced deviations are apparent (Left). 
Of course as t/T — > 1 these deviations would become in- 
creasingly pronounced also on the logarithmic scale. To 
asses the importance of correction terms for given values 
of a and T it is instructive to consider the ratio of (S 2 ) > 
for finite-time and infinite-time trajectories, 



R=(Sx 2 ) I (bx 2 ^ 



(7) 



(4) 



as function of the relative time-lag if = t/T. The re- 
sults for various cases are shown in Fig. |4] (Left). Thus, 
when a = 1.2 for instance, R(<p) decreases to 0.2 as if 
approaches 1. This finite-time depression is important 
in relating the amplitude of the measured time averaged 
MSD to the anomalous diffusion coefficient of the pro- 
cess. In the Brownian limit a — 2, R(ip) is independent 
of the finite measurement time T and equals 1. 

Finally we investigate the nature of the ergodicity 
breaking and, in particular, the role of the finiteness 
of trajectories. As already noted by Zumofen and 
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(a) 



3- 




(b) 




2.5 5.0 7.5 
r[10 5 A] 

FIG. 4: (a) R(tp) for various values of a and T; (b) Ergodicity 
breaking parameter as function of lag time, £B(t) for a = 1.2 
(circles) and a = 1.5 (squares) for T = 10 7 A (black) and 
T = 10 8 A (blue). The dashed and dotted red lines denote 
the theoretical values for infinite-time trajectories, and the 
full lines correspond to the predictions for the finite-time case. 



Klafter [§J the time averaged MSD differs from the cor- 
responding ensemble average, thus ergodicity is broken. 
We define the ergodicity-breaking parameter SB{t) = 
(fa 2 ) / (Ax 2 (t)) as ratio of time versus ensemble aver- 
aged MSD [31|. For our choice of tp(x,t) the ensem- 
ble averaged MSD asymptotically is (Ax 2 (r)) ~ 2(a — 
l)r 3 ~ a /((3-cv)(2-cv)) for 1< a < 2 as r -4 oo 6]. Thus 
£B(t) — I /(a — 1) as r — > oo, that is, time and ensemble 
averages differ only in terms of a constant. We call this 
phenomenon ultraweak ergodicity breaking, in contrast to 
the stronger weak ergodicity breaking of scale-free subd- 
iffusive processes. According to Eq. ([5]) we expect that 
the finiteness of trajectories will also affect £B(t). Inter- 
estingly £B(t) appears to be almost independent of r, as 
can bee seen in Fig. reffg4(b), but the value deviates sig- 
nificantly from 1/(q — 1) (dotted and dashed red lines). 
In fact, this is not surprising if considering the results 
in Fig. [3j where we found that the scaling of finite-time 
averages on the logarithmic scale agrees rather well with 
the prediction for infinite trajectories, suggesting that 
the correction terms effectively cause a rescaling of the 
generalized time-average diffusion coefficient. We might 
call this apparent ultraweak ergodicity breaking. 

We investigated the ultraweakly ergodic behavior of 
superdiffusive LWs, finding a pronounced scatter of ap- 
parent scaling exponents of the time averaged MSD for 
finite-time trajectories. These apparent scaling expo- 
nents range between ballistic motion (sticking to one 
velocity mode) down to subdiffusive values (localization 



due to erratic hopping between different velocity modes) . 
Moreover, averaged over many individual trajectories, 
the time averaged MSD is pronouncedly smaller than for 
very long trajectories. We quantify these effects in terms 
of an ergodicity breaking parameter. 

The present results reveal the importance to take into 
account the effects of the finiteness of trajectories when 
interpreting experimental results. They also demonstrate 
how the measured time series of different lengths re- 
veal more reliable information about the fundamental un- 
derlying dynamical process. The additional information 
comes from the dependence of time averaged quantities 
on the length of time series. Instead of attempting to 
measure or generate longer and longer time series to ex- 
tract reliable time averaged quantities, one could instead 
use many shorter time series and obtain even more re- 
liable results. Our results may also provide an alterna- 
tive and more robust method of determining exponents 
of probability densities of step durations. Namely, since 
we inevitably expect poor sampling of very long events 
this might be reflected in the obtained exponent. Using 
the time averaged MSD from measurements of different 
(but known) durations one should in principle be able to 
determine the exponent more accurately. 
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